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Abstract 

In quantum gravity, fields may lose quantum coherence by scattering off vac- 
uum fluctuations in which virtual black hole pairs appear and disappear. Although 
it is not possible to properly compute the scattering off such fluctuations, we argue 
that one can get useful qualitative results, which provide a guide to the possible 
effects of such scattering, by considering a quantum field on the C metric, which 
has the same topology as a virtual black hole pair. We study a scalar field on the 
Lorentzian C metric background, with the scalar field in the analytically-continued 
Euclidean vacuum state. We find that there are a finite number of particles at 
infinity in this state, contrary to recent claims made by Yi. Thus, this state is 
not determined by data at infinity, and there is loss of quantum coherence in this 
semi-classical calculation. 
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1 Introduction 



The possible loss of quantum coherence is one of the most exciting topics in quantum 
gravity. Recent work on D-branes has encouraged those that believe that the evaporation 
of black holes is a unitary process without loss of quantum coherence. It has been shown 
that collections of strings attached to D-branes with the same mass and gauge charges 
as nearly extreme black holes have a number of internal states that is the same function 
of the mass and gauge charges as e A ^ G , where A is the area of the horizon of the black 
hole [1, 2, 3]. They also seem to radiate various types of scalar particles [4, 5] at the same 
rate as the corresponding black holes. However, the D-brane calculations are valid only 
for weak coupling, at which string loops can be neglected. But at these weak couplings, 
the D-branes are definitely not black holes: there are no horizons, and the topology of 
spacetime is that of flat space. One can foliate such a spacetime with a family of non- 
intersecting surfaces of constant time. One can then evolve forward in time with the 
Hamiltonian and get a unitary transformation from the initial state to the final state. 
A unitary transformation would be a one to one mapping from the initial Hilbert space 
to the final Hilbert space. This would imply that there was no loss of information or 
quantum coherence. 

To get something that corresponds to a black hole, one has to increase the string 
coupling constant until it becomes strong. This means that string loops can no longer 
be neglected. However, it is argued that for gauge charges that correspond to extreme, 
or near extreme black holes, the number of internal states will be protected by non- 
renormalization theorems, and will remain the same. It is argued that there's no sign 
of a discontinuity as one increases the coupling, and therefore that the evolution should 
remain unitary. However, there's a very definite discontinuity when event horizons form: 
the Euclidean topology of spacetime will change from that of flat space, to something 
non-trivial. The change in topology will mean that any vector field that agrees with 
time translations at infinity, will necessarily have zeroes in the interior of the spacetime. 
In turn, this will mean that one cannot foliate spacetime with a family of time surfaces. 
If one tries, the surfaces will intersect at the zeroes of the vector field. One therefore 
cannot use the Hamiltonian to get a unitary evolution from an initial state to a final 
state. But if the evolution is not unitary, there will be loss of quantum coherence. An 
initial state that is a pure quantum state can evolve to a quantum state that is mixed. 
Another way of saying this is that the superscattering operator that maps initial density 
matrices to final density matrices will not factorise into the product of an S matrix and 
its adjoint [6]. This will happen because the zeroes of the time translation vector field 
indicate that there will be horizons in the Lorentzian section. Quantum states on such 
a background are not completely determined by their asymptotic behavior, which is the 
necessary and sufficient condition for the superscattering operator to factorise. 

One cannot just ignore topology and pretend one is in flat space. The recent progress 
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in duality in gravitational theories is based on non-trivial topology. One considers small 
perturbations about different vacuums of the product form M 4 x B, and shows that one 
gets equivalent Kaluza-Klein theories. But if one can have small perturbations about 
product metrics, one should also consider larger fluctuations that change the topology 
from the product form. Indeed, such non-product topologies are necessary to describe 
pair creation or annihilation of solitons like black holes or p-branes. 

It is often claimed that supergravity is just a low energy approximation to the fun- 
damental theory, which is string theory. However, the recent work on duality seems to 
be telling us that string theory, p-branes and supergravity are all on a similar footing. 
None of them is the whole picture; instead, they are valid in different, but overlapping, 
regions. There may be some fundamental theory from which they can all be derived as 
different approximations. Or it may be that theoretical physics is like a manifold that 
can't be covered by a single coordinate patch. Instead, we may have to use a collection 
of apparently different theories that are valid in different regions, but which agree on 
the overlaps. After all, we know from Goedel's theorem that even arithmetic can't be 
reduced to a single set of axioms. Why should theoretical physics be different? 

Even if there is a single formulation of the underlying fundamental theory, we don't 
have it yet. What is called string theory has a good loop expansion, but it is only 
perturbation theory about some background, generally flat space, so it will break down 
when the fluctuations become large enough to change the topology. Supergravity, on 
the other hand, is better at dealing with topological fluctuations, but it will probably 
diverge at some high number of loops. Such divergences don't mean that supergravity 
predicts infinite answers. It is just that it cannot predict beyond a certain degree of 
accuracy. But in that, it is no different from perturbative string theory. The string 
loop perturbation series almost certainly does not converge, but is only an asymptotic 
expansion. This means that higher order loop corrections get smaller at first. But after 
a certain order, the loop corrections will get bigger again. Thus at finite coupling, the 
string perturbation series will have only limited accuracy. 

We shall take the above as justification for discussing loss of quantum coherence in 
terms of general relativity or supergravity, rather than D-branes and strings. One might 
expect that loss of quantum coherence could occur not only in the evaporation of macro- 
scopic black holes, but on a microscopic level as well, because of topological fluctuations 
in the metric that can be interpreted as closed loops of virtual black holes [7]. Particles 
could fall into these virtual black holes, which would then radiate other particles. The 
emitted particles would be in a mixed quantum state because the presence of the black 
hole horizons will mean that a quantum state will not be determined completely by its 
behavior at infinity. It is with such loss of coherence through scattering off virtual black 
holes that this paper is concerned. Our primary intention is not to provide a rigorous 
demonstration that quantum coherence is lost, but rather to explore the effects that will 
arise, assuming that the semi-classical calculations are accurate, and it is lost. 



3 



In d dimensions, a single black hole has a Euclidean section with topology S d ~ 2 x R 2 . 
As has been seen in studies of black hole pair creation, a real or virtual loop of black 
holes has Euclidean topology S d ~ 2 x S 2 — {point}, where the point has been sent to 
infinity by a conformal transformation. For simplicity, we shall consider d — 4, but the 
treatment for higher d would be similar. 

On the manifold S 2 x S 2 — {point} one should consider Euclidean metrics that are 
asymptotic to flat space at infinity. Such metrics can be interpreted as closed loops of 
virtual black holes. Because they are off shell, they need not satisfy any field equations. 
They will contribute to the path integral, just as off shell loops of particles contribute to 
the path integral and produce measurable effects. The effect that we shall be concerned 
with for virtual black holes is loss of quantum coherence. This is a distinctive feature of 
such topological fluctuations that distinguishes them from ordinary unitary scattering, 
which is produced by fluctuations that do not change the topology. 

One can calculate scattering in an asymptotically Euclidean metric on S 2 x S 2 - 
{point}. One then weights with exp(— I) and integrates over all asymptotically Eu- 
clidean metrics. This would give the full scattering with all quantum corrections. How- 
ever, one can neither calculate the scattering in a general metric, nor integrate over all 
metrics. Instead, what we shall do in the next two sections is point out some qualitative 
features of the scattering in general metrics, that indicate that quantum coherence is lost. 
We shall then illustrate the effects of loss of quantum coherence and obtain an estimate 
of their magnitude by calculating the scattering in a specific metric on S 2 x S 2 — {point}, 
the C metric. It is sufficient to show that quantum coherence is lost in some metrics in 
the path integral, because the integral over other metrics cannot restore the quantum 
coherence lost in our examples. 

2 Lorentzian section 

We don't have much intuition for the behavior of Euclidean Green functions or their 
effect on scattering. However, if the Euclidean metric has a hypersurface orthogonal 
killing vector, it can be analytically continued to a real Lorentzian metric, in which it 
is much easier to see what is happening. We shall therefore consider scattering in such 
metrics. 

The Lorentzian section of an asymptotically Euclidean metric which has topology 
S 2 x S 2 — {point} will contain a pair of black holes that accelerate away from each 
other and go off to infinity. One might think that this is not very physical, but it is no 
different from a closed loop of a particle like an electron. Closed particle loops are really 
defined in Euclidean space. If one analytically continues them to Minkowski space, one 
gets a particle anti-particle pair accelerating away from each other. Any topologically 
non-trivial asymptotically Euclidean metric will appear to have solitons accelerating to 
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infinity in the Lorentzian section, but this does not mean that there are actual black 
holes at infinity, any more than there are runaway electrons and positrons with a virtual 
electron loop. One can regard the use of the Lorentzian metric, with its black holes 
accelerating to infinity, as just a mathematical trick to evaluate the scattering on the 
Euclidean solution. 

To understand the structure of these accelerating black hole metrics, it is helpful to 
draw Penrose diagrams. Start with the Penrose diagram for Rindler space with the left 
and right acceleration horizons, H al and H ar , and past and future null infinity, X~ and 
X + (see Figure 1). A uniformly accelerated particle moves on a world line that goes out 
to X~ and X + at the points where they intersect the acceleration horizons. One now 
replaces the accelerating particle and the similar accelerating particle on the other side 
with black holes. Thus, one replaces the regions of Rindler space to the right and left 
of the accelerating world lines with intersecting black hole horizons. It turns out that 
the two accelerating black holes are just the two sides of the same three dimensional 
wormhole, so one has to identify the two sides of the Penrose diagram, and the Penrose 
diagram will look like the one in Figure 2. At first sight it looks as if one has lost half 
of X~ and X + , but that is because this Penrose diagram applies only on the axis. One 
can get a better idea of the causal structure near infinity from Figure 3, in which a 
conformal transformation has been used to make X + into a cylinder S 2 x R 1 , with the 
null generators lying in the R 1 direction. The hypersurface orthogonal Killing vector 
of the Euclidean metric that allows continuation to a Lorentzian metric will be a boost 
Killing vector in the accelerating black hole metric and it will have two fixed points q 
and r onI + , lying on generators A and A' respectively. The past light cones of q and 
r minus the generators A and A' form the acceleration horizons. Thus one can see that 
nearly every null geodesic outside the black hole horizons goes out to X + in the region to 
the future of both acceleration horizons. The exceptions are the null geodesies that are 
exactly in the boost direction, which intersect the generators A and A'. We shall ignore 
A and A' as a set of measure zero on X + , and a number of the statements we shall make 
will be valid modulo this set of measure zero. 

3 Quantum state 

The analytically continued Euclidean Green functions will define a vacuum state \0)e 
which is the analogue of the so-called Hartle Hawking state [8] for a static black hole. 
The Euclidean quantum state can be characterized by saying that positive frequency 
means positive frequency with respect to the affine parameters on the horizons. In the 
accelerating black hole metrics there are two kinds of horizons, black hole and acceler- 
ation. Each kind of horizon consists of two intersecting null hypersurfaces, which we 
shall refer to as left and right, as in Figure 2. In choosing a Cauchy surface for the 
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spacetime (modulo a set of measure zero), we break the symmetry between left and 
right, and choose say the left acceleration horizon and the right black hole horizon. The 
quantum state defined by positive frequency with respect to the afline parameters on 
these horizons is the same as the quantum state defined by the other choice of horizons. 

Another Cauchy surface in the future (again modulo a set of measure zero) is formed 
by X + and the future parts of the black hole horizons H£ and H^, as in Figure 4. There 
is a natural notion of positive frequency on X + . On the black hole horizons the concept 
of positive frequency is less well defined. One could use Rindler time, but in any case, 
what one observes on X + is independent of the choice of positive frequency on the black 
hole horizons. 

The quantum state of a field on this background metric will be determined by data 
on either of these Cauchy surfaces. This means that the Hilbert space Ti. of quantum 
fields on this background metric will be isomorphic to the tensor products of the Fock 
spaces on their components: 



The vacuum state defined by the Euclidean Green functions is the product of the vacuum 
states of the Fock spaces for the left acceleration horizon and right black hole horizon; 



However, because of frequency mixing, the Euclidean quantum state won't be the prod- 
uct of the Fock vacuum states on J + and the future black hole horizons. Rather it will 
be a state containing pairs of particles. Both members of the pair may go out to X + , or 
both may fall into the holes, or one go out to X + and one fall in. 

Equation (1) shows that quantum field theory on an accelerating black hole back- 
ground does not satisfy the asymptotic completeness condition that the Hilbert space 
of the quantum fields on the background is isomorphic to the asymptotic Hilbert space 
of states on X + . Asymptotic completeness is the necessary and sufficient condition for 
scattering of quantum fields on the background to be unitary [6]. Thus there will be 
loss of quantum coherence. What happens is that to calculate the probability of observ- 
ing particles at X + , one has to trace out over all possibilities on the future black hole 
horizons. This reduces the Euclidean quantum state to what appears to be a mixed 
quantum state described by a density matrix. 

In a recent pair of papers [9, 10], Yi argued that the Euclidean quantum state in 
the Ernst metric would contain no radiation at infinity. The Ernst metric is similar to 
the metrics we are considering. However, in the explicit calculation that we carry out 
in the C metric, we find that there is indeed radiation at infinity. What's wrong with 
Yi's argument? As he was working with the Ernst metric, which isn't asymptotically 



= Fi+ ® T H + <g> T H + . 



(1) 



|o> B = |o) H JiW 



(2) 
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flat, he wasn't able to study the radiation at infinity directly. He therefore assumed 
that if there was no radiation on the acceleration horizon, there would be no radiation 
at infinity. But if we evolve some state forward from one of the acceleration horizons to 
Z+, part of the state can fall into the future black hole horizon. Therefore, there can be 
a non-trivial Bogoliubov transformation between the acceleration horizon and infinity, 
and Yi's assumption is incorrect. 

The Euclidean quantum state \0)e will be time symmetric, and so will contain both 
incoming and outgoing radiation. Unlike the Euclidean state for static black holes, there 
won't be radiation to infinity at a steady rate for an infinite time. Instead, the radiation 
will be peaked around the points q and r where the acceleration horizons intersect Z + . 
The radiation will die off at early and late times and the total energy radiated will be 
finite. 

Is this the appropriate quantum state? In the case of a static black hole, one usually 
imposes the boundary condition that there is no incoming radiation on Z~. This means 
that one has to subtract the incoming radiation from the Euclidean state to give what is 
called the Unruh state. This is singular on the past horizon, but that doesn't matter, as 
one normally replaces this region of the metric with the metric of a collapsing body. The 
energy for the steady rate of outgoing radiation comes from a slow decrease of the mass 
of the black hole formed by the collapse. However, in the case of a virtual black hole 
loop, there is no collapse process to remove the singularities on the past horizons of the 
black holes or supply the energy of the outgoing radiation. Therefore, we should study 
the Euclidean vacuum state, in which the energy of the outgoing radiation is supplied 
by the incoming radiation on I~. 

Our view therefore is that integrating over gauge equivalent virtual black hole metrics 
will cause the amplitude to be zero unless the energy of the outgoing particle or particles 
is matched by particles with the same energy falling in. One might object that one 
would never have exactly the combination of incoming particles that corresponded to 
the quantum state obtained from the Euclidean green functions. However, the Euclidean 
quantum state will appear to be a mixed quantum state on Z~ which contains every 
possible combination of incoming particles. One can choose one of these combinations 
as an initial pure quantum state that is incident on the virtual black hole loop. The final 
quantum state will then be that part of the Euclidean quantum state on Z + that has the 
same energy, momentum and angular momentum as the incoming state. Because of the 
trace over the future black hole horizon states, the final state on Z + will be mixed. Such 
an evolution from pure to mixed states can be described by a superscattering operator 
$ rather than an S matrix [6]. 

The dominant contribution will presumably come from virtual black hole loops of 
Planck size. The cross section for a low energy particle to fall into a Planck size static 
black hole is very low unless the particle is spin or 1/2 [11]. In the case of spin 1/2, the 
probability of emission will be reduced because the Fermi-Dirac factor (exp(cj/T) + 1) _1 
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tends to 1 at low uo while (exp(a>/T) — tends to T/uo. This suggests the effects 
of virtual black holes will be small except for scalar particles. In this paper we shall 
therefore do a scattering calculation for scalar particles in the C metric. This doesn't 
really qualify as a virtual black hole metric, because it has conical singularities on the 
axis, although one can interpret these as cosmic strings. We study the C metric because 
it has the same topological structure as a virtual black hole pair, but it has the great 
advantage that one can calculate the scattering, because the wave equation separates. 



4 C metric 



The charged C metric solution is [12] 

ds 2 = A~ 2 (x - y)~ 2 [G(y)dt 2 - G-\y)dy 2 + G~ 1 (x)dx 2 + G(x)d<p 2 

where 



(3) 



G(0 = (1 + r_AO(l - e - r + Af) = -r + r_A 2 (£ - - - &)(£ - ( 4 ) 



The gauge potential is 



(5) 



where q 2 = r + r_. We define m = (r+ + r_)/2. We constrain the parameters so that 
G(£) has four roots, which we label by £i < £2 < £3 < £4- To obtain the right signature, 
we restrict x to £3 < x < £4, and y to —00 < y < x. The inner black hole horizon lies at 
y — £1, the outer black hole horizon at y — £2, and the acceleration horizon at y — £3. 
The axis x = £4 points towards the other black hole, and the axis x = £3 points towards 
infinity. Spatial infinity is at x = y = £3, null and timelike infinity at x = y 7^ £3. This 
metric describes a pair of oppositely-charged black holes accelerating away from each 
other, although the coordinate system used in (3) only covers the neighborhood of one 
of the black holes. 

To avoid having a conical singularity between the two black holes, we choose 



47T 



(6) 



This implies that there will be a conical deficit along x = £3, with deficit angle 



S = 2?r 1 



(7) 



Physically, we imagine that this represents a cosmic string of mass per unit length 
(j, = S/8n along x = £3. At large spatial distances, that is, as x, y — > £3, the C metric 
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(3) reduces to flat space with conical deficit 5 in accelerated coordinates. The C metric 
also reduces to flat space if we set r + = r_ = 0. It reduces to a single static black hole 
if we set A = [13]. The limit r+A <C 1 is referred to as the point-particle limit, as in 
this limit the black hole is small on the scale set by the acceleration. 

The C metric was shown to be asymptotically flat in [14]. This is a considerable 
advantage over, say, the Ernst metric, as it means we will have a well-defined notion of 
X, and we can study the radiation at infinity directly. If we neglect the axis x = £3, all 
observers will intersect the acceleration horizon before reaching infinity, and the causal 
structure of the solution is roughly speaking given by the Penrose diagram shown in 
Figure 2. However, the metric is not spherically symmetric, so this diagram is not a true 
picture of the whole spacetime. We will refer to the left and right acceleration horizons 
as H a i and H ar , and to the left and right outer black hole horizons as and Hj, r . 
Further, the future and past halves of each horizon will be denoted by superscripts ±. 
Hopefully the diagram clarifies the meaning of this notation. 

We will only discuss the behavior at future null infinity. As the metric is time- 
symmetric, the discussion of past null infinity will be identical. We can conformally 
compactify the C metric by using a conformal factor Q = A(x — y). The conformally 
rescaled metric is 

ds = Q 2 ds 2 = G{y)dt 2 - G-\y)dy 2 + G'\x)dx 2 + G{x)d V 2 . (8) 

Null infinity is the surface Q = 0, that is, x = y (more precisely, its maximal extension; 
the coordinate system of (8) misses the generator on which the other black hole intersects 
X + [14]). The induced metric on X + is 

dsl = G{y){dt 2 + d<p 2 ). (9) 

Note that, at null infinity, t is a spatial coordinate. The normal to X + is 

n a = V a n = 2AG(y)d y . (10) 

We see that t and tp are constant along the orbits of n a , which are the generators of X + , 
so they are good coordinates on the manifold of orbits of X + . It is convenient to define 
new coordinates 9, 77 where 

(so Arj = 2ir). We also make a further conformal rescaling with a conformal factor 
Q' = \G'fa)\ sm9/2G^ 2 (y), so that 

dsl = n' 2 dsl = d9 2 + sin 2 9dif '. (12) 
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Figure 2: A Penrose diagram for the C metric, neglecting the axis x = £3. The heavy 
dashed lines are singularities, and the surfaces are boundaries of the spacetime. A 
Cauchy surface C for the region outside the inner black hole horizons constructed from 
one black hole horizon and one acceleration horizon is shown. 
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Figure 3: The structure of X + in the conformal gauge (12). The two points are where the 
black holes intersect X + , and their past light cones are the acceleration horizons. Two 
of the 9, r) cross-sections are pictured. The dashed lines represent the conical deficits in 
the metric (3); they are not part of X + . 

In this conformal gauge, an affine parameter along the generators of X + is 



which labels the 9,r] cross-sections. The structure of X + in the conformal gauge (12) 
is depicted in Figure 3. In this conformal gauge, X + is divergence-free, and 9,7] are 
coordinates on the manifold of generators of X + , so we can see that X + has topology 



We can obtain the Euclidean section of the C metric by setting t = ir in (3). To make 
the Euclidean metric positive definite, we need to restrict the range of y to £2 < V < £3- 




(13) 



It is also useful to define another coordinate 




(14) 



S 2 x R. 
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There are then potentially conical singularities at y = 6 and y = £3, which have to 
be eliminated. We can avoid having a conical singularity at y = £ 3 by taking r to be 
periodic with period 

AT = " = ^- (16) 

In this paper, we assume the black holes are non-extreme, that is, 6 < 6- We can then 
only avoid having a conical singularity at y = £2 by taking the two horizons to have the 
same temperature, so that both conical singularities can be removed by the same choice 
of At. This implies 

6-6 = 6 -£ 3 . (16) 

The Euclidean section has topology S 2 x S 2 — {pt}. This Euclidean section can be used 
to study the pair creation of black holes by breaking cosmic strings [15, 16, 17]. However, 
we want to use it simply to determine the appropriate vacuum state on the Lorentzian 
section. Since the black hole and acceleration horizon have the same temperature on 
the Euclidean section, the analytic continuation will give Green's functions which are 
thermal with temperature 1/(5 with respect to the time parameter t in the Lorentzian 
section. 

The region of the spacetime outside the inner horizon of the black holes is globally 
hyperbolic. Consider a Cauchy surface for this region which is made up of one black 
hole horizon and one acceleration horizon (say the left acceleration horizon and the right 
black hole horizon), as pictured in Figure 2. As explained earlier, the Hilbert space is 
isomorphic to the tensor product of the Fock spaces on the two horizons (1). Positive 
frequency on the Fock spaces is defined with respect to the affine parameter along the 
horizon. The state we wish to study is the analytically-continued Euclidean vacuum 
state \0)e given in (2). 

In the next section, we will describe the solution of the scalar wave equation on the 
C metric background. We then use this to calculate the Bogoliubov coefficients in the 
subsequent section. 



5 Scalar Wave Equation 

We consider a minimally-coupled massless neutral scalar field, so the wave equation 
is just D0 = 0. One of the great advantages of considering the C metric is that this 
equation separates. It is easy to see this if we observe that the C metric is a solution of the 
vacuum Einstein-Maxwell equations, and hence R = 0. The minimally coupled equation 
above is therefore equivalent to the conformally-invariant equation □</> — = 0. 
But in solving this latter equation, we are free to make conformal transformations. In 
particular, we can transform to the conformal gauge (8), in which this equation takes 
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the form 

(17) 

where because of the conformal rescaling, = <f>/A(x — y). Thus we see that if we use 
the ansatz 

<P = A(x-y)e^ t e^v(x) 1 (y), (18) 
then we get two second-order ODEs for v{x) and 7(1/), 

d x [G(x)d x u(x)} - -^rn 2 u(x) + [\d 2 x G{x) + D]v{x) = (19) 

and 

d y [G{y)d yl {y)\ + ^fHv) + [\d 2 y G(y) + £>] 7 (y) = 0, (20) 

where D is a separation constant, and G(£) is given in (4). Note that tp is a periodic 
coordinate with period 47r/|G"(£ 4 )|. Thus m = m \G'(^4)\/2, where m is an integer. 
We assume, without loss of generality, that it is positive. 

One way to rewrite these equations that offers some further insight is to define new 
coordinates 

z= Iwr x =lwr (21) 

which have the advantage that d z = G(y)d y ,d x = G(x)d x . Note that the integral for 
z in (21) diverges as we approach a horizon, as G(y) — * at the horizons. Thus, 
— 00 < z < 00 only covers the region between two of the horizons; similarly, £3 < x < £4 
is mapped to —00 < x < 00 ■ We can write (19,20) as 

d 2 x u(x( X )) - m 2 u(x( X )) + V eff (xHx( X )) = 0, (22) 

d 2 zl {y{z)) + cA(y(*)) + V eff {zMy(z)) = 0. (23) 

That is, (20) reduces to the one-dimensional wave equation with effective potential 
V e ff(z), which is given by 

V eff (z) = G(y(z))^d 2 G(y(z)) + D]. (24) 

There is a similar expression for V e ff(x)- It is not possible to invert (21) to obtain y(z) 
explicitly, but we can make some observations. Near the horizons, G(y) — > 0, and thus 
the effective potential becomes unimportant, so 7(2/) ~ e ±luJz . Similarly, near x = £3,^4, 
v(x) ~ e ±mx . Obviously, for physically-interesting solutions, we must have v(x) ~ e~ m ' x l 
as x ~~ * ±00. 
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We can rewrite the metric (3) in terms of these coordinates: 

ds 2 = A~\x - y y 2 [G(y)(dt 2 - dz 2 ) + G(x)(d X 2 + dip 2 )], (25) 

where by x,y we mean x(x), y( z )- This coordinate system evidently only covers the 
region between two of the horizons (or between the acceleration horizon and infinity). 
That is, there is a coordinate system like this for each of the diamond-shaped regions 
in the Penrose diagram in Figure 2. We will therefore refer to these as the Rindlerian 
coordinates. We can now define null coordinates u, v — t ± z. Since z increases as we go 
from the acceleration horizon towards the black hole horizon, the u and v coordinates 
run as shown in Figure 4. Thus, u is a (non-afline) parameter along H^ r and Hj^, while v 
is a (non-afline) parameter along and H^. As is usual for bifurcate Killing horizons, 
these parameters are related to the affine parameters U, V on the acceleration horizon 
by u — - In \U\, v — —\ In | V|, where k = G"(£ 3 )/2 is the common surface gravity of the 
two horizons. 

These coordinates are useful for specifying boundary conditions near the black hole 
and the acceleration horizons, and we will see later that we can easily write down explicit 
forms for the positive-frequency wavefunctions on the horizons in terms of them. How- 
ever, as we can't write V e ff explicitly as a function of z, we can't solve the differential 
equations in this form. 

If we return to the initial forms (19,20) for the ODEs, we find that they can be 
considerably simplified. In the simplification, we will exploit the equal-temperature 
condition (16), which imposes an additional symmetry on the form of G(£). If we make 
a coordinate transformation 

*-te^)K-5fc+«], (26) 

then 

G(0 = -|(?-« 2 )(?-l), (27) 

where 

C = r+r.x4-6)" Q = t^i- *=5<*-6>. < 28 > 
Note that a > 1, £, ip > 0, and that = ipd^. If y and x are defined in terms of y and 
x following (26), then the inner black hole horizon is at y = —a, the outer black hole 
horizon is at y — — 1, and the acceleration horizon is at y — 1, while the range of x is 
1 < x < a. In terms of these coordinates, 

5fG(0 = ~[12f -2(l + a 2 )], (29) 
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Figure 4: A Cauchy surface C for the region outside the inner black hole horizons 
constructed from X + and the future halves of the black hole horizons. The Rindlerian 
coordinates u, v between the acceleration and outer black hole horizons are also shown. 



16 



so it is convenient to define 



1 

6 ( 



/3f> = -(l + a*H 



2 ' 



so that 



^G(0 + J> = 



We can now write z explicitly; 

dy 



-J 



c 



G{y) 2(a 2 - 1) 





a + y 


+ In 


i-y 




a 


a-y 




l + y 





(30) 
(31) 

(32) 



We can now see clearly that z diverges at the event horizons y = — a, ±1. We can further 
see that z — > — oo as we approach y = —a, 1, the inner black hole and acceleration 
horizons, and z — > oo as we approach y = — 1, the outer black hole horizon. There is a 
similar explicit expression for x, and x ^ -oo as we approach x — 1 and % — > oo as we 
approach x — a. The consideration of the form (22,23) suggests a further simplifying 
transformation. If we set 



u(x) = e mx h{x) 



a + x \ 2a (° M ) 



x 



a — x 



and 



7(y) 



7(y) 



1 + x 



n(£) 



2(c 2 -l) „ 



then we can finally rewrite (19,20) as 



d £ [(x 2 - l){x 2 - a 2 )d £ n{x)] - 2m(dzh{x) + 2{x 2 - (3 2 D )n{x) 



l){y 2 - a 2 )dyf{y)} + 2iu<^/(y) + 2(y 2 - (3 2 D )f{y) 



0, 
0. 



(33) 
(34) 

(35) 
(36) 



This is the simplest form in which we can write these equations. 

We have been able to simplify the form of the wave equation considerably. However, 
(35,36) still have five regular singular points, at £ = ±1, ±a, oo, so they can't be solved 
exactly. We will therefore need to use some further simplifying assumption in solving 
the wave equation. There is only one dimensionless parameter in the metric, r + A, 
as the equal-temperature condition fixes r^A as a function of r + A. Therefore we are 
driven to consider the point-particle limit r + A <C 1. In this limit, a ~ 1 + 4r + A, and 
( ~ 8r + A PS 2 (a — 1). For reasons of convenience, we will use (a — 1) as the small 
parameter. 
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6 Bogoliubov Transformations 



Having laid the groundwork, we can now define and evaluate the Bogoliubov coefficients. 
We can write the field operator in terms of annihilation and creation operators on the 
Hilbert spaces associated with the black hole and acceleration horizons: 

<P = A(X- y)^ m / MfLti* m + RlJ&m + ZlmKlm + fLfiln), W 

where /* Jm , f^ im are sets of positive frequency modes which have non-zero support on the 
black hole and acceleration horizons respectively, b^ lm , b® lm are the particle annihilation 
operators, and b b j lm , 6°| m are the particle creation operators. Here, positive frequency 
means with respect to the affine parameters U, V on the horizons. 

Following [18], we see that a suitable set of positive frequency states on the black 
hole horizon is 

tiim = |1 _ e -L /K|1/2 ^ m (^)[^ + e-™^], (38) 

where v\ m is a solution of (19) with D given by (3d = 1 + 21(1 + 1), and are functions 
which are non-zero on the future and past parts of the black hole horizon respectively, and 
which are positive frequency with respect to the Rindler parameter, that is, = e~ luu . 
We know already that only a discrete set of values for m are allowed, and we will see 
below that the same is true for /. We wish to normalize the modes so that (f^i m , f b >i> m ') = 
5 mm '5u'5(uj — w'), which implies |iV| 2 = l/(47r|u;| A</?). Note that although the positive- 
frequency solutions are labeled by a frequency uj, they do not have a single frequency 
with respect to U, and the solutions are still wholly positive frequency with respect to 
U when uj is negative. For this to be a complete set of positive-frequency solutions, we 
must allow uj to run over — oo < uj < oo. One can write down a similar set of positive 
frequency solutions on the acceleration horizon. 

In appendix A, we consider (35) with (a — 1) <^ 1, and we learn that, as we might 
have expected, there is a restriction on the form of the data on the black hole horizon. If 
we write I — lo + 0(a — 1), then the solutions v\ m (x) will only be regular at both of the 
axes x = £3, x = £4 if Iq is an integer and Iq > m , where mo is the integer appearing in m. 
In the point-particle limit, the x, tp section approaches spherical symmetry, so Iq is the 
usual total angular momentum quantum number, while m is the angular momentum 
with respect to the axis along which the black holes are accelerating. 

We can also write the field operator in terms of modes which are positive frequency 
at infinity: 

= A(x -y) J duj d9 dr] (p^ + p^al + + q u cl), (39) 

where p u are a set of modes with non-zero support on X + which are positive frequency 
with respect to f, and a u , are the corresponding annihilation and creation operators. 
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The modes have non-zero support on the future black hole horizon, and c u , ct are 
the corresponding annihilation and creation operators. We won't bother to define these 
latter modes, as their form is irrelevant to the calculation of particle production on X + . 
Following [19], we choose the positive frequency modes to have the form 

Pu = 7^=z . fl 6(6 - e )5( v - r) Q ) = ' &» 6(e - e )S( v - Vo ) (40) 

V 2ttcj sin \>2-kAAuj 

on X + , in the conformal gauge where the metric on X + has the form (12). We define 
oj = (I>|G"(£ 4 )| sin6> /4A Each mode is thus non-zero on one generator of X + , labeled by 
#0,7/0, and has frequency oj with respect to the afline parameter along that generator. 
The complete set of positive frequency modes is given by < u < oo. They are 
normalized so that (pw,p' w ) = 2u;<5 3 (A; — k'), where k is the three-momentum, and points 
in the direction (9q,t]q). 

Since both sets of modes are complete bases for the space of solutions of the wave 
equation, we can write one in terms of the other. That is, 

ft'lm = J duj d9 dr) (a^imPw + fiLsimPu + terms involving q u ), (41) 

and similarly for j^, Xm . If we substitute these expansion into (37), and require consistency 
with (39), then we find that 

O-ui = J dw'( a lu/'lmPl'lrn + filui'lmPj'lm + a L' Irrfiu 1 Im + fiwu/'lmpj'lm)- (42) 

The quantities a h ^, lw /3^ m , cOj ro , fil^'im are called the Bogoliubov coefficients. Since 
we know how the annihilation and creation operators which were defined on the horizons 
act on |0) Ei to determine how the annihilation and creation operators defined at infinity 
act on \0)e, we just need to compute these coefficients. 

The operator we are most interested in is the number operator N u = a^a^, which 
gives the number of particles in the mode p w . In the state \0)e, 

(0\N U \0) E = X lM J Wdcy / (^ l ^, I , m ,(0|^ I X'iw|0>s + ...) (43) 

/oo 
du'd Pics' lm\ 2 + \fiuui'lm\ 2 )i 
-oo 

where we have expanded by (42), and in the second line we have used the canonical 
commutation relations and the fact that b^, lm \0)E = 0,b°, lm \0)E = 0. 

We should now calculate the Bogoliubov coefficients fil^'im and fiZu'im- However, it 
turns out to be quite difficult to calculate the latter coefficient. Therefore, we wish to 
argue that it is sufficient to calculate the contribution from the Bogoliubov coefficient 
associated with the black hole horizon filvn m ; the other contribution should be similar. 
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We broke the symmetry between the left and right horizons when we defined the 
Euclidean vacuum state, by defining it to be the product of the vacua of the Fock spaces 
for the left acceleration horizon and right black hole horizon. However, the vacuum state 
is in fact symmetric under left-right interchange. That is, it is also equal to the product 
of the vacua of the Fock spaces for the right acceleration horizon and left black hole 
horizon. Take the vacuum state on C and evolve it forward through the right Rindler 
diamond, from H~ t and H br to H+. and H bl . There will then be correlations between 
H br and if+, due to the correlations between the two halves of the black hole horizon in 
the Cauchy surface C. Further, there are no correlations between H br and because 
on C, the state has no correlations between the black hole and acceleration horizons. 
Since the state is left-right symmetric, the correlations between the two halves of the 
acceleration horizon in the Cauchy surface C can therefore only give rise to correlations 
between H bl and H^, and these correlations will be related to the ones coming from the 
black hole horizon. Both sets of correlations give rise to correlations between J + and 
the future black hole horizons, which give the particle creation, so the particle creation 
due to the acceleration horizon should just be the image under the left-right interchange 
of the particle creation due to the black hole horizon. This justifies our only calculating 
the latter contribution. 

We now calculate fl^u'im- The modes v^iV^i^A^ are orthogonal, and 

(p w ,jv) = 2u)8 3 (k - k') = ^—-5^ - r]' )S(9o - 9' )5(u> - Cj 1 ). (44) 

to sin uq 

Thus, we can use (41) to show 

- b _ ^sinflp b 2Au b 

Pwui'lm — 2 Ju'lm) — ^Q/^^\PoJi Jaj'lm)- ) 

To evaluate this inner product, we need to express both the modes as functions on the 
same Cauchy surface. We do this by evolving the mode f^i m forwards from C to C. 

The propagation from C to C can be broken up into two stages: propagation through 
the right Rindler diamond, from and H b ~ to H+. and H bl , and propagation through 
the future diamond, from and if+ to J + . The initial data on H br is just the 
restriction of (38) to the past part of the black hole horizon, while f^, lm vanishes on H~ t . 
From the discussion of (23), we recall that at the acceleration and black hole horizons, 
7(2/) ~ e ±tuJz . Using this and the form (38) of the mode f^i m , we find that the boundary 
conditions on 7(1/) are 

1 (y)= e - wz + C R e wz (46) 
at the black hole horizon z — > 00, and 

7(1/) = C T e^ z (47) 
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at the acceleration horizon z — > — oo, where Cr and Ct are constants which remain to 
be determined. In appendix B, we solve (36) with these boundary conditions in the limit 
r + A <C 1, assuming oo ~ 0(1), and find that Ct ~ (a — l) 2l+1 , and that, for l = 0, 
\Ct\ ~ (a — l)u/2. Because the transmission factor Ct is increasingly suppressed for 
increasing /, we will be mostly interested in the contribution from the Iq — mode, as the 
other contributions will be smaller than the terms that we neglect in our approximate 
calculation of the Iq — contribution. 

The propagation from and if + to X + is also described in appendix B. This part 
of the calculation is substantially easier; it is very similar to solving the angular equation 
(35). In the conformal frame where the metric has the form (12), the restriction of f^ im 
to X + is 

fb I , — 2NC T G 1/2 (y) -icj'(t+z)im<p-m\z\f (\~ (\ (ao\ 

Ju,'lm\l+ ~ |j _ e -2nw>/K\l/2\Q,^\ sin e 6 6 Jlu>'\P) n lm{P)- 

In this expression, /j<y(p) is given by the definition at the end of appendix B, and we 
have defined nj m (p) to be ni m (p) for p < 1/2 (x < 0) and e 2mx ni m (p) for p > 1/2 (x > 0), 
where ni m (p) is the approximate solution of the angular equation given in appendix A. 
When p — > 0, fi UJ i{p)- l fii m {p) — > 1. When p — > 1, fii m {p) — > e lg , some constant phase. 
Evaluating the inner product, we find that 



PL'im = ^- —— / dze iu,r e iu, ' z e- m ^f lu>/ (p)n lm (p), (49) 



where t — t corresponds to 6 — 9 , if — ip corresponds to i] = i] , and we have used 
dr = dz/G 1 ^ 2 (y), which follows from (14) and (21). 

Note that apart from an overall phase, this expression depends only on the frequency 
u>, and not on 9 , r] . This means that the expression is boost invariant, that is, invariant 
under translations in t, as the orbits of the boosts are the cross-sections labeled by r, 
and thus these boosts preserve the frequency uo with respect to r. 

We can't evaluate the integral in (49), but we can still get some interesting physical 
information about the radiation out of this expression. Because G(y) — > as z — > ±oo, 

dv ^ 

- = - ±oo when z -+ ±oo, (50) 

and hence the e lujr part of the integrand oscillates with an effective frequency which 
tends to infinity at large \z\. Since the amplitude is bounded, the main contribution to 
the integrand will come from the region near z = where the integrand oscillates slowly. 

The integral in (49) will give an answer which is peaked in uj' with some finite width, 
so the integration over u/ of |/3£ w // m | 2 in (43) should give a finite answer. By contrast, 
in the case of a static black hole, the analogous formula for the Bogoliubov coefficient 
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gives a delta function in uj', so the expected number of particles is infinite (that is, in 
that case there is a steady flux of particles across X + ). 

Our calculation of the transmission factor in appendix B is only valid for \uj'\ < 1, 
and we might expect that for sufficiently large uj', the potential barrier would become 
unimportant, and Ct ~ 0(1). However, the Bogoliubov coefficient will be small for 
large negative uj' because of the factor |1 — e~ 2lvu) / K \~ 1 / 2 . We also expect that it would 
be small at large positive uj', as the integrand in the integral in (49) will then oscillate 
rapidly for all values of z, making the integral small. Thus, the main contribution to the 
integral over uj' in (43) will come from small negative uj' , where the calculation of Ct is 
valid. 

We expect that the size of the contribution from each l,m will be primarily deter- 
mined by the transmission factor, so we expect that the contribution from l = m = 
will dominate the summation over /, m in (43). We now consider the form of this contri- 
bution in the point-particle limit, where we can somewhat simplify the expressions and 
illustrate some of these remarks. When (a — 1) <C 1, we have G(y) ~ 4p(l — p) on X + , 
where p — (y — l)/(a — 1). Further, z ~ \ ln(p/ (1 — p)), as < p < 1 on X + , 



so 



G (y) « — rs~- ( 51 ) 

cosh z 

Thus, dr/dz ps coshz, and hence 

r ss sinhz. (52) 
Also, k ~ 1, fooj(p) ~ 1, and n o(p) ~ 1- Therefore 

PL'oo « - A ; ; / ^ e *(^-^*). (53) 

As we argued above, the main contribution to the integration will come from the region 
near z = 0, so the primary contribution to P^'ooi an d hence to the number operator, 
will come from the part of the generator closest to the points where the black holes 
intersect X + . If we restrict our attention to the region near z = 0, we can expand sinhz 
in a power series, and we see that the integrand is most nearly constant near z — if 
io' = —cu, so we expect that the Bogoliubov coefficient will be peaked at uj' = —uj. This 
peak will become narrower as uj — > 0, approaching a delta function in the limit, but the 
amplitude tends to zero in this limit because of the factor of uj in front of the integral, 
so this does not imply infinite particle production. 

The leading-order part of the total particle production along the generator labeled 
by 0o,r} is given by integrating |/5^/ 00 | 2 over uj and uj'; we can't do this integral, but 
given the arguments above, it seems reasonable to expect the answer to be finite. The 
integration over all generators, which gives the total particle production, will not give 
rise to any divergences either. 
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7 Discussion 



In the first part of this paper, we argued that the scattering off virtual black hole 
pairs, which could lead to loss of quantum coherence in ordinary scattering processes, 
could be discussed in terms of a path integral over Euclidean metrics with topology 
S 2 x S 2 — {point}. In this approach, one considers the scattering in each metric and 
performs a path integral over all such metrics. Since we cannot perform this path 
integral, we then restricted the discussion to one such metric, and analytically continued 
the solution to a Lorentzian section to make the scattering easier to understand. 

We argued that the appropriate quantum state is the analytically-continued Eu- 
clidean vacuum state |0)_g, and we argued that this state will contain a finite, non-zero 
number of particles at infinity. It is well-known that from the point of view of an observer 
co-moving with the black holes, this state corresponds to a thermal equilibrium between 
the black holes and a thermal bath of acceleration radiation. Thus, this state must be 
time-reversal invariant, which means that the particle content at past null infinity X~ is 
the time-reverse of the particle content at future null infinity X + . This implies that no 
net energy is gained or lost by the black holes in this scattering process, which is what 
we would expect for a model of a virtual loop, and is in agreement with the fact that 
the state is an equilibrium as seen by co-moving observers. 

The fact that there is a non-zero number of particles at X + implies that there is loss of 
quantum coherence in this semi-classical calculation, as each particle detected at infinity 
can be thought of as one member of a virtual pair, the other one of which has fallen 
into the black hole, carrying away information. More formally, there are correlations 
between modes on future infinity and modes on the future black hole horizon, and the 
information encoded in these correlations is lost because we do not observe the state on 
the future black hole horizon. This loss of quantum coherence is of the same character 
as that observed in static black holes. 

In the second part of the paper, we proceeded to an explicit calculation of the 
scattering in the C metric. Although the Euclidean C metric solution has topology 
S 2 x S 2 — {point}, it is not usually thought of as describing a virtual black hole loop, 
as it is a solution of the field equations, and it has a conical singularity along one of the 
axes. However, we believe it is a reasonably good model for a virtual black hole loop, 
and the wave equation separates in this background, so it is relatively easy to study the 
scattering explicitly. The C metric is asymptotically flat [14], so it is also straightfor- 
ward to study the radiation at infinity. One slightly surprising fact about the structure 
at infinity is that the affine parameter along generators of X + is f, which is spacelike 
between the black hole and acceleration horizons, while the boost time coordinate t 
becomes a spacelike coordinate labeling the generators of T + . 

It is also worth noting that the transmission factor Ct ~ (a — l) 2l+1 . This implies 
that the dominant contribution to the particle production is in the s-wave, as for static 
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black holes, because of the high centrifugal potential barrier for higher-spin modes. It 
also suggests that the scattering of higher-spin fields off such virtual black hole loops 
will be suppressed relative to that of scalar fields, as they cannot radiate in the s-wave. 
This is in agreement with the arguments of [20, 19]. 

The calculation we have actually been able to perform is rather limited; we considered 
only one specific, rather special metric, and we were only able to study the scattering 
on it in a particular limit. However, the results we have obtained give an estimate of 
the magnitude and nature of the effects of virtual black hole loops, and they agree well 
with our general expectations. 
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A The angular quantization condition 

In the point-particle limit (a - 1) < 1, the deviations from spherical symmetry in the 
x, ip part of the metric become small, so we would expect that the dependence on x will 
reduce to the usual angular momentum modes, with quantum numbers / and m. Recall 
that because of the periodicity of if, m = m |G'(^4)|/2 = m [l + 0(a — 1)], where m 

is an integer. We also expand I — l + l±(a — 1) + The range of x is 1 < x < a, so 

we define a new coordinate p = (x — l)/(a — 1). If we expand n(x) in powers of a — 1, 
n(x) = riimip) — n o(p) + (« — l)ni(p) + . . ., then (35) can be separated into a series of 
equations for these functions. The first equation is 

d P [p{p - l)d p n (p)} - m d p n (p) - l (l + l)n (p) = 0. (54) 

This equation is a hypergeometric equation. The possible values of Zo are restricted 
by requiring that the solution behave appropriately at the two poles, p — 0,1. As we 
said earlier, for the solution for to be physically relevant, we must have v(x) ~ e~ m ' x ' 
as x ^ ±oo. That is, we require that v{x), and hence 0, doesn't blow up at the 
axes. Since v(x) = e mx n(x), the appropriate boundary conditions on ni m (p) are that 
nim(p) = 1 as x — oo, which corresponds to p — 0, and rii m (p) ~ e _2mx as x ^ oo, 
which corresponds to p — 1. Therefore, the appropriate solution of the hypergeometric 
equation (54) is no(p) = F(l + 1, — 1 ; 1 + m ;p), where F is the hypergeometric series, 
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as F(a, b; c; p) — > 1 as p — > 0. If we analytically continue this solution to a neighborhood 
of p — 1, we find 

/ \ r(l + m )r(m ) i \ 

w) = w TvF7 1 / ,^ F(l + l,-lo,l-rno',l-P) 

T(m - l )T(m + l + l) 

+ r rr/ + r?lw~M ) ( 1 - ^) m ° F ( m o ~ M + ™o + /o; 1 + m ; 1 - p). (55) 

1 (/ + ijl { — Lq) 

The second term has the appropriate behavior for p — > 1, since e~ x ~ (1 — p) 1 / 2 for 
p « 1. Thus, the coefficient of the first term must vanish, which can only happen 
if Iq — m is a non-negative integer. This is just the usual quantisation condition for 
angular momentum, and Iq is thus the total angular momentum quantum number. 

The next-order term l\ can similarly be fixed by requiring that the solution n\{p) is 
regular at p — 0, 1. Unfortunately, it is not possible to give a general formula for l\\ the 
equation must be solved separately for each Iq, itlq. We are particularly interested in the 
case Iq = m = 0, as we expect this mode to make the dominant contribution to the 
particle production on X + . In this case, uq{p) = F(l, 0; l;p) = 1, while the equation for 
ni(p) is 

d p \p(p - l)9 p m(p)] =k-p. (56) 

This equation has a solution which is regular at p — 0, 1 only if l\ = 1/2; in this case, 
the solution is ni{p) = —p/2 + C , where C is a constant. One can similarly fix all the U. 



B The transmission factor 

In section 6, we found that to evolve the positive-frequency modes from C to C, we need 
to calculate the transmission factor Ct between the black hole and acceleration horizons. 
That is, we need to solve (36) with the boundary conditions (46,47), and find Ct- For 
convenience, we will repeat those here. The equation is 

dy[(f - l)(f - a 2 )dyf(y)] + 2iuCdyf{y) + 2(f - (3 2 D )f(y) = 0. (57) 

In terms of the function f(y), the boundary conditions are 

f(y) = l + C R e 2 ^ (58) 

near the black hole horizon y — — 1 and 

f(y) = C T (59) 

near the acceleration horizon y — 1. 
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We can't solve this equation exactly, but if (a — 1) <C 1, then we can solve it 
approximately. First note that if y 2 — 1 is 0(1) (that is, if y is not close to ±1), we can 
neglect terms involving a — 1 to approximate (57) as 

dylif ~ l)%fm + 2(y 2 - Pl)f(y) = 0. (60) 

In neglecting the term involving u, we have made the further assumption that \u\ ~ 0(1); 
that is, that oo is not large. This equation is now a hypergeometric equation. To put 
it in the standard form, we set f(y) = 2 a (l — y 2 )~ a (a — l) a g(s), where s = (y + l)/2, 
a — I + 1. Then 

s(s - l)d 2 s g(s) - 2/(2s - l)d 9 g(s) + 2/(2/ + l)g(s) = 0, (61) 

where we have used (3d = l + 2/(/ + l). We use / rather than Z in the approximate 
equations in this section, because regarding I as an integer would introduce degeneracies 
in the approximate equations which are not present in the exact equation. Near y — ±1, 
the solutions of (61) can be expressed in terms of hypergeometric series about y = ±1. 
However, we cannot approximate (57) by (61) in a neighborhood of radius 0(a — 1) 
around y — ±1, which is precisely where we wish to impose boundary conditions. 

Therefore we need a separate approximation to cover these neighborhoods. When 
y 2 — 1 ~ (a — 1), make a coordinate transformation y = ±(1 + (a — 1)?±)- Then if we 
keep just the leading terms, (57) becomes 

d q± [q±(q± - l)d q± f(q±)] ± iwd q J - 1(1 + 1)/ = 0, (62) 

where f(q±) = f(y)- These are, once again, hypergeometric equations. The solution 
about y — — 1 which satisfies the boundary condition (58) is 

f(q_) = F(a, b;2-c; q_) + C R (-q^F(b + c - 1, a + c - 1; c; (63) 

and the solution about y — 1 which satisfies the boundary condition (59) is 

f(q + ) = C T F(a,b;c;q + ), (64) 

where F is the hypergeometric function, a — I + 1, b — — I and c = 1 — iuo. Now 
analytically extend these solutions to large q±: at large q_, the solution (63) becomes 

, T(c)T(a - b) ( Cr + r(2 - c)r(c - b)\ , _ ,_ b 



T(c-b)T(a) \ n T(c)T(2-c-b) 
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while at large q + , the solution (64) becomes 

Ha \ -C r(c)r(5 ~ a) r „ )-« + C T ^ a - b h „ Y b (66) 

Now for 1 <C 1 1 <C (a — l) -1 , both approximations are applicable, so we can use 
the large-distance behavior (65,66) of the approximation for y near ±1 as boundary data 
for the approximation (61). If we pick the solution g(s) to be 

a(s) ~ r ( c ) r ( fe ~ a ) (r | T(2-c)T(c-a) \ 

9{s) ~ r( c - a)r(b) \ Cr + r(c)r(2 -c-a)) F{ ~ 21 ' ~ 21 " ~ 2 '' s) (67) 
, r(c)r( fl - 6) ^ + r(2 - c )r(c-6) \ (q _ lf -a T - bs a- bm 1; 2 , + 2; a)> 



r( c -&)r(a) v r(c)r(2- c -6) 

then the boundary conditions obtained from (65) are automatically satisfied. We can 
analytically continue this solution to a neighborhood of s = 1; to satisfy the boundary 
conditions obtained from (66) in this neighborhood at the same time, we must require 

T(c)T(b-a) _ T(c)T(a-b) ( r(2 - c)T(c - b) \ b b 

Ct T(c - a)T(b) ~ r(c - b)T(a) [ Cr + r(c)r(2 - c - b) ) ( " " 1} 2 (68) 

and 

r(c)r(a-6) 6 _ _ 6 _ r(c)r(6-a) / r(2- c )r( c - a ) \ 

CT r(c-6)r(a) (a_1) 2 -T(c-a)T(b) [ Cr + T(c)T(2 - c - a) ) ' (69) 
Solving these two equations for Cr and Ct, we find 

^ = -e'^ (70) 

and 

C R = -e^ + 5C T , (71) 

where 



„ T(2-c)T(c-b) 

r(c)r(2- c -6) 

and 

'a-l\ B - 6 r(&-a)r(a)r(c-6) 



r(a-6)r(6)r(c-a)' 
Note that these coefficients satisfy |Ct| 2 + |C/?| 2 = 1, as they should. 



(72) 
(73) 
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After some manipulation, we find 



4i /a-l\ 2l+1 T(l + l-iu)T(l + l + iuj 



5-5 = - — ; - ( — — ) — W7 - M2 -sinh7ra;. (74) 

2> 



21 + 1 V 8 ) r(z + 



Also, 5 ~ (a — l) 2l+l , so the denominator in Ct can be ignored for this leading-order 
calculation. For large /, we thus find 



C T « 2eW+?U ^j^j sinhTTcu, (75) 
while for l = 0, we find 

Ct « c «*-5) f — ) w . (76) 



2 

These results are valid for (a — 1) <C 1 and |u;| < 1. 

We have found the value of f(y) at the acceleration horizon y = 1. The region 
between if^, and X + is the region between y = 1 and y = x; to evolve f(y) through 
this region, we just need to find the form of f(y) between y = 1 and y = a, which will 
also be the solution on X + . Now, the approximation (64) is valid throughout this region, 
so the result is simply that on X + , 

f(y) = C T fiM « C T F(a,b;c;p), (77) 

where a, b, c are as in (64). Note that x = y implies q + = p. For Iq = 0, the leading-order 
part of this solution is /o w (p) ~ 1, just as for ni m (p). 
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